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Abstract. In this second part we prove that the full nonlinear fluid-solid 
system introduced in Part I is stabilizable by deformations of the solid that 
have to satisfy nonlinear constraints. Some of these constraints arc physical 
and guarantee the self-propelled nature of the solid. The proof is based on the 
boundary feedback stabilization of the linearized system. Prom this boundary 
feedback operator we construct a deformation of the solid which satisfies the 
aforementioned constraints and which stabilizes the non linear system. The 
proof is made by a fixed point method. 

1. Introduction. In this introduction, let us begin with completing and presenting 
again the fluid-solid model which has been introduced in the first part of this work, 
so the reader can skip this presentation. 

1.1. Presentation of the model. The movement of the solid in the inertial frame 
of reference is described through the time by a Lagrangian mapping denoted by Xg, 
so we have 

S{t) ^ Xs{S{0),t), t>0. 

The mapping Xs{-,t) can be decomposed as follows 

Xs{y,t) = hit) + Ilit)X*{y,t), ye 5(0), 

where the vector h{t) describes the position of the center of mass and R(t) is 
the rotation associated with the angular velocity of the solid. In dimension 3 the 
angular velocity is a vector field whereas it is only a scalar function in dimension 

2. However, can be immersed in R"^ and this scalar function can be read on 
the third component of a 3D-vector. More generally in this work, since all the 
calculations made in dimension 3 make sense in dimension 2, we will consider only 
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vector fields of M.^ and matrix fields of M'^'^'^. For instance, w and R are related to 
each other through the following Cauchy problem 

— - bHR ^ with§H = 







UJ2 


r 





-UJl 









W2 = 0. 







I R(0) = li 

where in dimension 2 we have uj = uj^ and uJi 
The couple {h{t), R(t)) describes the position of the solid and is unknown, whereas 
the mapping X*{-,t) can be imposed. This latter represents the deformation of the 
solid in its own frame of reference and is considered as the control function on which 
we can act physically. When this Lagrangian mapping X*{-,t) is invertible, we can 
link to it an Eulerian velocity w* through the following Cauchy problem 

dX* 

-^{y,t)^w*{X*{y,t),t), X*{y,0) = y-h{0), yeS{0). 

If Y*{-,t) denotes the inverse of X*{-,t), we have 

f)X* 

w*{x*,t) = __(r*(x*,t),i), X* es*it)^x*is{0),t). 

The fluid flow is described by its velocity u and its pressure p which are assumed 
to satisfy the incompressible Navier-Stokes equations. For X* satisfying a set of 
hypotheses given further, the system which governs the dynamics between the fluid 
and the solid is the following 

du 

— -i/Au + (it- v)m + Vp = o, xeJ'it), te(o,oo), (i) 
ot 

divu^O, xeJ^it), <e(0,oo), (2) 

w = 0, xedO, te(0,cx3), (3) 

u=h'(t)+Lu{t)A{x~h{t))+w{x,t), xedS{t), te(0,oo), (4) 



Mh"{t)^- a{u,p)ndT, te(0,oo), (5) 

Jas{t) 

(luj)' (<) = -/ {x- h{t)) A a{u,p)ndT, t £ (0, oo), (6) 

JdS(t) 

u(y,0) = ?/o(y), y e ^(0), /i'(0) = /n e R^^ w(0)=woeR^ (7) 



where 



S{t) ^h{t)+Ii{t)X*{S{Q),t), Tit) ^0\S{t), (8) 

and where the velocity w is defined by the following change of frame 

w{x,t) = R(t) (R(t)^(a;-/i(<)),t) , a; G (9) 

Without loss of generality, we assume that h{0) = 0, for a sake of simplicity. The 
symbol A denotes the cross product in R^. The linear map lo A • can be represented 
by the matrix S{uj). In equations (5) and (6), the mass of the solid M is constant, 
whereas the inertia moment depends a priori on time. In dimension 2 the inertia 
moment is a scalar function which can be read on the inertia matrix given by 

i{t) = n^^psix,t)\x-h{t)\'dxji^3. 
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In dimension 3 it is a tensor written as 

I{t) = [ ps{x,t){\x-h{t)\H^s-{x^h{t))®{x^h{t)))dx. 
Js{t} 

The quantity ps denotes the density of the soUd, and obeys the principle of mass 
conservation 

where VXg denotes the Jacobian matrix of the mapping ■ For a sake of simphcity 
we assume that the sohd is homogeneous at time t = 0: 

Ps{y,0) = ps>0- 

In system (l)-(9), v is the kinematic viscosity of the fluid and the normalized 
vector n is the normal at dS{t) exterior to It is a coupled system between the 

incompressible Navier-Stokes equations and the differential equations (5)-(6) given 
by the Newton's laws. The coupling is in particular made in the fluid-structure 
interface, through the equality of velocities (4) and through the Cauchy stress tensor 
given by 

cr(w, p) = 2vD{u) -pld = v {Vu + {Vuf^ - p Id. 

Indeed, the Dirichlet condition (4) partially imposed by the deformation of the solid 
(through the velocity w) influences the behavior of the fluid whose the response is 
the quantity a{u,p)n in the fluid-solid interface. It represents the force that the 
fluid applies on the solid, and then it determines the global dynamics of the solid 
(through equations (5) et (6)) and thus its position. 

The problem is the following: What is the deformation X* of the solid we have to 
impose in order to stabilize the environing fluid and thus induce a behavior of the 
fluid which stabilizes the velocities of the solid? We shall assume a set of hypotheses 
on the control function X* , that we state as follows: 

HI: For ah t > 0, X*{-,t) is a C^-diffeomorphism from 5(0) onto S*{t). 

H2: In order to respect the incompressibility condition given by (2), the volume 
of the whole solid has to be preserved through the time. That is equivalent 
to say that 

/ — — • (cofVX*) ndF ^ 0. (10) 

H3: The linear momentum of the solid is preserved through the time, which 
leads to 

P5(y,0)X*(y,0dy = 0. (11) 

5(0) 

H4: The angular momentum of the solid is preserved through the time, which 
leads to 

BX* 

ps{y,Q)X*{y,t)A—-{y,t)Ay = 0. (12) 

5(0) Ot 
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Imposing constraints (17) and (18) enables us to get the two following constraints 
on the undulatory velocity w 



As equations (5) and (6) are written, the constraints (13) and (14) are implicitly 
satisfied in system (l)-(9). Hypotheses H3 and H4 are made to guarantee the 
self-propelled nature of the motion of the solid, that means no other help than its 
own deformation enables it to interact and to move in the surrounding fluid. 

1.2. The main result and the methods. The existence of global strong solutions 
for system (l)-(9) has been proven in [18] in dimension 2, and more recently in 
dimension 3 for small data in a framework where the deformations of the solid are 
limited in regularity. The main result of this second part is Theorem 7.1, that we 
can state as follows: 

Theorem 1.1. Let be A > 0. Then for {uq, ft-i, cjq) small enough in 

there exists a deformation X* satisfying the hypotheses H1-H4 given above and 

also 



The proof of this theorem is based on the preliminary stabilization of the lin- 
earized system in Part I. The idea is the following: If the perturbations of the system 
(which are represented by the initial conditions) are small enough, the behavior of 
the nonlinear system is close to the evolution of the linearized system. The strat- 
egy we follow in this second part consists in rewriting system (l)-(9) in cylindrical 
domains (that means in domains whose the space component does not depend on 
time), by defining a change of variables and a change of unknowns. Then we focus 
on the nonlinear system thus obtained. The feedback boundary control obtained 
in Part I enables us to stabilize the nonhomogeneous linear part of this system, 
whereas we have now to consider as a control function a deformation of the solid 
which satisfies the nonlinear constraints H1-H4 given above. So the difficulty is 
to define properly from this boundary feedback control a deformation of the solid 
which stabilizes the full nonlinear system. For that the boundary feedback control 
is extended inside the solid as it is done in the section 7 of Part I; The deformation 
thus obtained satisfies the linearized version of the constraints stated in the hy- 
potheses H1~H4. Then we project the displacement associated with this mapping 
on a set representing the displacements satisfying the nonlinear constraints. The 
deformation obtained is said to be admissible^ that is to say it lies in a well-chosen 
functional space and it obeys the hypotheses H1-H4 (see Definition 2.1). This pro- 
jection method enables us to decompose the deformation velocity - on the fluid-solid 
interface - into two parts: The first part of this decomposition stabilizes the linear 
part of the nonlinear system, whereas the remaining part satisfies good Lipschitz 
properties with respect to the boundary feedback control. This point is essential if 
we want to prove by a fixed point method that such a deformation stabilizes the 




(13) 




ps{x,t){x - h{t)) hw[x,t)dy = 0. 



(14) 




< C. 
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nonlinear system. 

Besides the technical aspects of this work, a particular difficulty is the consideration 
of a control that has to satisfy nonlinear physical constraint. The originality of our 
contribution can be read in a perspective which concerns the study and the control 
of the swim of a deformable structure at an intermediate Reynolds number. By 
controlling the velocity of the environing fluid in a bounded domain, the solid sta- 
bilizes to zero the full system which is already dissipative (because of the viscosity), 
but the strength lies in the fact that this stabilization is obtained for all prescribed 
exponential decay rate. 

Other papers treat of this issue: Let us quote for instance the work of Khapalov 
[11, 12] where the incompressible Navier-Stokes equations are considered; More re- 
cently let us mention the work of Chambrion & Munnier [5, 6] dealing with perfect 
fluids, and where geometric methods have been used; Concerning the swim at a 
low Reynolds number, let us quote the recent paper [13] where prescribed types 
of deformations are considered. At a high Reynolds number, the work of Glass & 
Rosier [9] consists in applying the Coron return method in order to prove the local 
controllability of the position and the velocity of a boat, which is able to impose a 
velocity on a part of its boundary, in order to move itself in a fluid satisfying the 
incompressible Euler equations. The result we give concerns only the velocities of 
the system. 

1.3. Plan. The functional framework is given in section 2 for the unknowns as 
well as for the control function and the changes of variables. Then the change of 
variables and the change of unknowns are introduced in section 3. Technical results 
in relation with the changes of variables are stated and proven in Appendix A and 
Appendix B. In section 3 we also give the nonlinear system written in cylindrical 
domains. The main result of Part I is recalled in section 4 where the feedback 
stabilization of the nonhomogeneous linear system is treated. Then we construct in 
section 5 an admissible deformation which is supposed to stabilize the full nonlinear 
system in section 6. The link with the actual unknowns is made in section 7 which 
can be viewed as a conclusion. 

2. Definitions and notation. We denote by = -^(0) the domain occupied by 
the fluid at time t = 0, and by 5 = 5(0) the domain occupied by the solid at t = 0. 



We assume that S is smooth enough. We set S*{t) = X*{S,t), F*{t) =0\ S*{t), 
5^=5x(0,cx3), Q^=.Fx (0,oo), 



For more simplicity, we assume that the density ps at time < = is constant with 
respect to the space: 



Ps{y,0) = PS >0- 

We assume without loss of generality that ho = 0. This implies in particular 



and 



Qoo = U J"(t) X {t}. 



t>0 
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The cofactor matrix associated with some matrix field A is denoted by cof A. Let 
us keep in mind that when this matrix is invertible we have the property 

cofA^ = (detA)A-i. 

Let us introduce some functional spaces. We use the notation 

^ [R'{n)f or [ff (17)]^' 

for some integer k, for all bounded domain ft of or R'^. 

2.1. Functional framework for the unknowns. The velocity u will be consid- 
ered in the following functional space 

a^'^Qoo) = L2(0,oo;H2(^(i)))nHi(0,cx);L2(.F(i))) 

that we endow and define with the norm given by 

2 



du , . 



dt < +O0. 

L2(.F(t)) 



Analogously we can define spaces of type H^i (0, T; H"^ (f](t))) and (0, T; H"^ {^(t))) 
for all time-depending domain fi(i), where si and S2 arc non-negative integers. The 
pressure p will be considered in L^(0, T; H^(J^(i))); At each time t it is determined 
up to a constant that we fix such that Jj-^^^p = 0. Thus in particular from the 
Poincare-Wirtinger inequality the pressures P defined in T can be estimated in 
Hi (J") as follows^ 

\\P\\uH^} < C||VP||l^(^). 

The same estimate will be considered for other functions which play the role of a 
pressure in J-{0). 

More classically in the cylindrical domain Q'^ we set 

and we keep in mind the continuous embedding 

a^HQl) L-(0,oo;Hi(J-)). 

2.2. Functional framework for the control and the changes of variables. 

Let us introduce some functional spaces for the solid displacements - defined in the 
domain iS - and the changes of variables which will be defined in the domain J-. 
We will mainly consider mappings X* satisfying X* {■,()) = IdjF-, and thus we will 
consider the displacements 

Z*^X*-lds e Wa(^^), 
where the space 'W\{S^) is defined as follows 



fl7* 



dZ* 

Z*(2/,0)-0, — (y,0) = \fyeS. 



^In the following the symbols C or Co will denote more generally a generic positive constant 
which does not depend on time or on the unknowns. 
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We endow it with the scalar product 
(^1 ; ^2)w^(SS„) 



dt 



dt 



L2(0,oo;H3(5))nHi(0,cx3;Hi(5)) 

which makes it be a Hilbert space, because of the continuous embedding 

Wa(^L) L°°(0,oo;H3(5))nWi-°°(0,oo;Hi(5)). 
Indeed, for X* — lAg G yVA('5'So)i have the foUowing estimates 



(15) 



dX* 



dt 



L~(0,oo;Hi(S)) 



< 



< 



< 



,-As 



dt 



ds. 



H3(5) 



dX* 



dt 



L2(0,oo;H3(5)) 



-As 



at2 



ds, 



Hi(5) 



1 + A 



At 



dX* 



dt 



Hi(0,oo;Hi(S)) 



In the rewriting of system (l)-(9) in cyhndrical domains, and in the final fixed 
point method, the mapping which has the most important role is denoted by X. 
We could consider the same type of functional space for this mapping in Qj^, but 
we will only need to consider - and we will only have - estimates oi X — Idjr in the 
space 

L°°(0, cx); H3(J-)) n Wi'°°(0, cx); Hi(J')). 

Thus for more clarity we set 

UiSl^) = l2(0,oo;H3(5)) nHi(0,c5o;Hi(5)), 
H{Ql,) = L2(0,oo;H3(J-))nHi(0,cx;Hi(J-)), 
VV(Q^) = L°°(0,oo;H3(J-))nWi'°°(0,oo;Hi(J-)). 

The main reason for which we choose to consider the aforementioned Lagrangian 
mappings in such functional spaces is the following: The changes of variables X will 
be - indirectly - obtained through extensions of the deformations of the solid X*] 
We will need to consider displacements X — lAjr which have the regularity indicated 
in the spaces given above, and for that we will have to consider displacements of 
the solid which lie in a Hilbert space (for the projection method) and which satisfy 
at least 

.xtdX* 



dt \ds 



e l2(0,oo;H5/2(55)) nHi(0,oo;Hi/2^55)). 



2.3. Definitions of stabilization for the nonlinear system. Let us specify the 
notion of admissible control and the definition we give to the stabilization of the 
main system. 

Definition 2.1. Let be A > 0. A deformation X* \s said admissible for the nonlinear 
system (l)-(9) if the displacement obeys 

X*-lAs e Wx{S°^), 
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if X*{-,t) is a C^-diffeomorphism from S onto S*{t) for all t >0, and if for alH > 
it satisfies the following hypotheses 

0, (16) 
0, (17) 

0. (18) 

Remark 1. First, the constraint which forces X* to be a C^-diffeomorphism can be 
relaxed if the control X* stays close to the identity Ids- Indeed, in this work the data 
arc assumed to be small enough, so that we will lead to consider the displacement 
X* - Ids small enough in the space yV\{S'^) and so in L°°(0, oo; H^(5)); Thus we 
can assume that this constraint is always satisfied. 

Definition 2.2. We say that system (l)-(9) is stabilizable with an arbitrary expo- 
nential decay rate if for all A > there exists an it admissible deformation X* (in 
the sense of Definition 2.1) and a positive constant C - depending only on uq, hi 
and Wo - such that the solution (u,p, h! of system (l)-(9) satisfies 

l|e'^*(u,P, ^','^)||H2.i(Qoo)xL2(0,oo;Hi(J^(t)))xHi(0,oo;R3)xHi(0,oo;R3) < C. 

3. Change of variables and change of unknowns. In order to use a change of 
unknowns which will enable us to rewrite the main system (l)-(9), we first extend 
to the whole domain O the mappings Xsi-^t)^ initially defined on S. Then by this 
extension - denoted hy X - we will define new unknowns whose interest lies in the 
fact that they are supposed to satisfy a system written in cylindrical domains, that 
is to say domains whose the space component does not depend on time. 

3.1. The change of variables. For a vector field h E 11^(0, oo; K^) and a rotation 
R G H^(0,oo;M^) which provides an angular velocity u! e H^(0, oo; R"^), and for an 
admissible deformation X* - in the sense of Definition 2.1 - we construct in this 
subsection a mapping X which has to satisfy the properties 

detVX = 1, inJ'x(0,oo), 
X = h + RX*, ona5x(0,oo), 
X = Iddo, on do x (0,oo), 

and such that for all t > the function X{-,t) maps onto J-{t), dS onto dS{t), 
and let invariant the boundary dO. The construction of such a mapping is quite 
technical, that is why we develop the details of this construction only in Appendix A, 
in the same time as the regularity deduced on X. However, let us note an important 
point: 

Remark 2. In Appendix A, the definition of the mapping X is conditioned by a 
smallness assumption on the solid deformation X* . Since the deformation which will 
be chosen in sections 6 and 7 for stabilizing the full nonlinear system will actually 
depend only on - and be controlled by - the initial data (uo,/ii,'^o) that we will 
assume small enough, it is possible to proceed like this. 

Note that the definition of such a change of variables is also done in [18] or in 
[7]. The change of variables we use in this paper has the same properties as the 
ones constructed in these articles. But the way we proceed here is not the same: 
In [18] the mapping X is constructed by extending the Eulerian velocity w to the 



—-iy,t)-{cofVX*iy,t))ndr{y) 
dS ot 

' X*{y,t)dy 



X*{y,t)A—^{y,t)dy 
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fluid part, but this means is not suitable in a framework where the role of the 
Lagrangian mapping representing the deformation of the solid is central and where 
its regularity is limited; Indeed the velocity w{-,t) is defined through w*{-,t) which 
is itself defined on the domain X*{S,t). Concerning the means used in [7], the 
problem solved for constructing the mapping X in this paper is similar to ours, but 
it requires a smallness assumption on the time existence. Here the hypothesis we 
make is the smallness of e"^* in an infinite time horizon space. 

3.2. Rewriting the main system in cylindrical domains. Let us transform 
system (l)-(9) into a system which deals with non-depending time domains. For 
that we set the change of unknowns 



uiy, t) = R(t)^u(X(2/, t),t), u{x, t) = R(i)u(r(x, t), t), 

Piy, t) = p{X{y, t),t), p{x, t) = p{Yix, t), t), 



(19) 



for X £ F{t) and y £ and 

h'{t) = R(i)'^/i'(t), Cj{t) = R(i)^w(t). (20) 

Remark 3. Let us notice that if h' and uj are given, then by using the second 
equality of (20) we see that R satisfies the Cauchy problem 

^(R) = §(Rw)R = R§(w) ■,^,^|~^ ( ° 
(\t with bi(ijj) — \ ujj, —UJ 

R(< = 0) = Iks, \—Cj2 uji 



So R is determined in a unique way. Thus it is obvious to see that, in (20), h' and 
w are also determined in a unique way. Moreover, since we have 

u{x, t) = Ii{t)u{Y{x, t),t), p{x, t) = p{Y{x, t),t), 

and since the mapping Y depends only on /i, a; and the control X*, we finally see 
that if (w,p, h',uj) is given, then {u,p,h' ,uj) is determined in a unique way. 

For what follows, it is convenient to define the mappings 

X{y, t) = R{t)^{Xiy, t) - hit)), {y, t) e T x (0, oo), 

I I ^ (21) 
Y{i,t) = Y{h{t)+R.{t)S:,t), {i,t)e\jR{t)^{T{t)-h{t))x{t}. 

t>a 

The regularity, dependance with respect to the unknowns, and estimates for the 
mappings X and Y are given in Appendix B. 

Then, like in [7], system (l)-(9) whose the unknowns are {u,p, h',uj) is rewritten 
in the cylindrical domain T x (0, oo) as the following system whose the unknowns 
are (w,p, h' ,uj): 

dii ~ 

— - v'Lu + M(u, /i',a;) + Nit + w(f) A it + Gp = 0, in J" x (0,oo), 

div u = g, in X (0, oo), 

it = 0, in do x (0, oo), 

dX* 

u = h'{t) + w{t) A X*{y,t) + ^{y,t), (y, t) € 95 x (0, oo). 
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Mh"{t) = d{u,P)\/Y{XYndT - MCo{t) A h'{t), t G (0, oo) 

JdS 

I*{t)Q'it) = - / X*iy,t) A (aiu,p)\/Y{Xfn) dV 
'95 ^ ^ 

-I*'{t)Cj{t) + r{t)Cj{t) A Cj{t), t e (0, oo) 



where specifies tlie i-tli component of a vector 
[LuUy,t) ^ [Vu{y,t)AY{Xiy,t),t)], + \/^u,{y,t) : (vYVY^) {X{y,t),t), 



Miu,h',u){y4) = -^u{y,t)^Y(X(y,t),t) i^' {t) + i::{t) A X{y,t) + ^{y.tyj 

^u{y,t) - Vu{y,t)^Y{X{y,t),t)u{y,t), 
Gpiy,t) =VY{X{y,t),tfVpiy,t), 

d{u,p){y, t)^v (Vu(y, t)yY{X{y, t),t) + VY{X{y, t), tfyu{y, tf) - p{y, t)h 

and 

g{y,t) = trace (Vw(y,i)(lR3-Vy(x(y,t),t) 

Notice tliat we liave 

cof(Vl) = det{WX)WY{Xf = WY{Xf, 

and so from the Piola identity we can actuaUy express this nonhomogeneous diver- 
gence term as 5 = div G, where 

Giy,t) - (lK3-Vy(X(y,t),t))u(y,t). 

For A > 0, we now set the following change of unknowns: 

u = e^*w, p = e^*p, h' = e^*h', Co ^ e^*Cj. (22) 

The idea of this second change of unknowns is the following: If we find a control X* 
such that the quadruplet {u,p, h' , uj) is bounded in some infinite-time horizon space, 
then the intermediate unknowns {u^p, h' ,uj) will be stabilized with an exponential 
decay rate, and it will be sufficient to deduce from that the same property for the 
actual unknowns (it,p, h',Lu) (see section 7). 

The system satisfied by (u,p, h',uj) is then transformed into 

— - Xu- i^Au + \/p^ F{u,p,h',u}), in J" X (0,00), (23) 

div u = div G{u), in J" x (0, 00), (24) 

u = 0, on X (0,00), (25) 

dX* 

u = h'{t)+uj{t)Ay + e^*^— + W{u}), (y, <) G (95 x (0, 00), (26) 

ot 
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Mh"-XMh' = -j (j{u,p)ndr + FM{u,p,h',uj), in(0,c»), (27) 

Jas 



IqCj' {t) - XIquj ^ - I y Aa{u,p)ndr + Fj{u,p,uj), in(0,oo), (28) 

JdS 

u{y,0) = uo{y),yeJ', /i'(0) = /ii G M^ w(0) = wo e M^ (29) 

with 

F{u,p,h',u}) = j^(L-A)M-e-^*M(M,/i',w)-e'^*Nu-(G-V)p-e"^*wAu, 
G{u,h\Lo) = (lR3-Vy(X(j/,i),t))^i, 
W{lj) = ljA{X*^ Id) , 
FM{u,P,h',(2]) = -Me-^^u Ah'it) 

-V j ^Vu (vf(X) - Irs) + {\/Y{X)^l^,^^ ViF^VY{XfnAT 

- f aiu,p)hY{X)-l^sY ndV, 
JdS ^ ' 

Fi{u,p,uj) = - {I* - Io)uj' + X{F* - Io)lv - r'cj + e-^^ruj Auj 

-V j yA {^u (vf (X) - Irs) + (Vy(X) - Ir3)^Vm^) Vf (l)^ndr 

yA (<j{u,p){VY{X)-l^.) 



OS 



dr 



+ / {X* - Id) A {a*{u,p)\/Y{Xfn^ dr. 



as 

Remark 4. An important remark is the fohowing: Since system (l)-(9) and the 
system satisfied by {u,p, h' ,u}) above are equivalent, and since for an admissible 
control X* satisfying the constraint (16) the compatibility condition is satisfied for 
system (l)-(9), in system (23)-(29) the underlying compatibility condition enables 
us to have automatically the following equality 

/ G{u,h',Lo)-ndT = / (e^*^^ + W{Lb)) -ndV 
Jas Jas \ ot J 

as soon as u = on dO. 



4. The stabilized nonhomogeneous linearized system. 

4.1. A nonhomogeneous linear system. In this section, let us consider the 
nonhomogeneous linear system suggested by the writing of system (23)-(29): 

— - AJ7 - i/AC/ + VP = F, inJ-x(0,oo), (30) 
ot 

div U = div G, in T x (0, oo), (31) 



C/ = 0, in aO X (0,oo), (32) 
U = H'it) + n{t) Ay + (:iy,t)+W{y,t), ina5x(0,oo), (33) 
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MH"{t) - XMH'{t) = - [ a{U, P)ndT + ¥m, t e (0, oo), (34) 

JdS 

lo^'it) - Xloilit) ^ - [ y Aa{U,P)ndr + ie(0,oo), (35) 

JdS 

U{y^O)^UQ{y), yeT, i?'(0) = /ij G R^ ^(0) = £ M^ (36) 
We assume that the data satisfy 
Fe L2(0,c3o;L2(J-))^ 

Ge H2.1(Q[^), We L2(0,c3o;H3/2(a5))nHi(0,oo;Hi/2(a5)), 

Fm e L2(0,c»;R3)^ F/ e L2(0,oo;R3), 

and the foUowing compatibihty conditions 

C-ndr = 0, [ G-ndT^ [ W-ndT, G(-,0)=0, = 0. 

OS JdS JdS 

Remark 5. Note that in this system the control C does not stand for e'^*^^, 
but it represents a boundary control which satisfies the linearized version of the 
constraints (16)-(18). That is why this notation C is only used for a control for a 
linear system, like in Part I. Thus the right-hand-side term W in (33) stands for 

4.2. Operator formulation. The nonhomogeneous divergence condition (31) can 
be lifted by setting 

U = U-G. 

Then the quadruplet {U , P, H', Cl) is supposed to satisfy the following system 

dU - . . 

— - AC/- J/A[/ + VP = F + Fg, in J- X (0,cx3), 
div C/ = 0, in J" X (0, od), 

C7 = 0, in5Ox(0,oo), 
U = H'{t) + n{t) Ay + (:{y,t) +W{y,t) -G{y,t), in55x(0,oo), (37) 

MH"{t) - \MH\t) = - [ cr(C7, P)ndT + Fm + Fm,g, t S (0, oo), 

JdS 

rn'it) - XloClit) = - [ yA cr(U, P)ndr + Fj + Fj.c, t € (0, oo), 
JdS 

U{y, 0) = uo(2/), yeJ^, iJ'(O) = /ii e M^ f2(0) = wo e M', 

where 

Fg - -?^ + AG + i.AG, 
ot 



Fm,g - -2i^ / DiG)ndT, 

JdS 



F/.G = -2iy / yAD{G)ndT. 

JdS 

By following the steps of the operator formulation used in Part I for the homoge- 
neous linear system, the right-hand-side W — G in (37) can be lifted and the pressure 
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of the resulting system can be eliminated, so that the latter can be formally rewrit- 
ten as follows 

^ = AjJ + BxC + F + Bx{W-G) (38) 
(Id-P)C7 = (Id-P) (io(^') + ^o(!^)) 

with U= (¥U,H\ny and F = (F + Fg, Fm + Fm.g, Fj + F/,g). In this formula- 
tion the operators P. Lq, A\ and B\ arc given in Part I (sec section 4). 

Let us now recall the most important result established in Part I of this work. 

4.3. Boundary feedback stabilization. We replace the control ( by the feedback 
operator denoted by IC\ and defined in section 6 of Part I. Then in the evolution 
equation (38) the operator A\ becomes A\ + B\IC\. This latter is stable, so that 
we can estimate 

II'^IIh2,1(Q!^) + ||-P||l2(0,oo;H1(J^)) + l|-ff'llHl(0,oo;K3) + 

< C* (IIuoIIhi(JF') + |^i|r3 + I^oIrs + i|P||L2(0,oo;L2(JP)) + ||G||h2.1(Q0^) 

+ l|W'||L2(0,oo;H3/2(JF))nHl(0,oo;Hl/2(jr)) + 1 1 Fm 1 1 L^ (0,oo;R3) + 1 1 1^/ 1 1 L^ (0,oo;K3 ) ) • 

Let us keep in mind that we have U = U + G, and so we have also the following 
estimate 

II^I1h2.i(Q!L) + ll-P|lL2(0,oo;Hi(Jf)) + l|-ff'llHi(0,oo:R3) + jlf^!lHi(0,oo;K3) 

< C (!|mo||hi(jc-) + \hi\m3 + IcjoIrs + ||P||l2(o,oo;L2(j^)) + 1I'G1|h2,i(qo^) 

+ l|W||L2(0,oo:H3/2(jr))nHi(0,oo:Hi/2(jr)) + 1 1 F^/ || l2 (o,oo;R3) + 1 1 IF/ || L2 (0,.oo;R3) ) • 

(39) 

for some independent constant C > 0. 

5. Construction of a stabilizing admissible deformation. Let us consider a 
boundary stabilizing control ( e L'^{0,oo;'H.^/^{S)) n ll^{0,oo;H^/^{S)) which can 
be chosen in a feedback form, as described above. The purpose of this section con- 
sists in defining from this boundary function a deformation X* which is admissible 
in the sense of Definition 2.1, and which has to have a satisfying Lipschitz behavior 
with respect to the function C. 

5.1. statement. The main result of this section is the following: 
Theorem 5.1. Let he Q e L2(0, oo; H5/2(5)) n Hi(0, cx); Hi/2(5)) satisfying 

I C-ndT = 0. 
Jas 

If C is small enough in L2(0, oo; H5/2(5)) n Hi(0, c»; Hi/2(5)), then th ere exists a 
mapping X* which is admissible in the sense of Definition 2.1, and which satisfies 

11-'^* - Id5||wA(SS„) < llCllL2(HV2(5))nHi(Hi/2(5)), (40) 
= O (||C||L2(H5/2(5))nHi(Hi/2(5))) . 

L2(H5/2(a5))nHi(Hi/2(a5)) 

(41) 



At 



dX* 
~dt 
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Moreover, if two functions Ci o-nd Q2 are close enough to in L^(0, 00; H^/^(5)) H 
H^(0, 00; H^/^(iS)), then the admissible deformations and X2 that they define 
respectively satisfy 

\\^2 ~ ^iWwxiS^) < IIC2 - Cl|lL2(H5/2(5))nHi(Hi/2(5)), (42) 



L2(H5/2(a5))nHi(Hi/2(a5)) 



< |jC2 - Cl|lL2(HS/2(5))nHi(Hi/2(5)) 

X-f'^ (IICl|lL2(HV2(5))nHi(Hi/2(5)) + |lC2||L2(H=/2(S))nHi(Hi/2(5))) , (43) 

with K{r) — > when r goes to 0. 

The proof of this theorem is divided into two parts. In the first one, from a 
boundary funetion C, we construct an internal solid deformation satisfying the 
linearized versions of constraints (16)-(18). Then in a second time we project 
the displacement X'^ — Ids on a space of displacements which define admissible 
deformations. 

5.2. First step of the proof. The first step of the proof of Theorem 5.1 consists 
in defining from a boundary control C an internal deformation X^ which satisfies 
the linearized versions of constraints (16)-(18). For that, let us recall a result of 
Part I, which is a consequence of the addition of Proposition 6 and Corollary 1 of 
section 7: 

Proposition 1. For C e L2(0, cx); H5/2(5)) n H1(0, cx); Hi/2(5)) satisfying 

C-nAT = 0, 



Ids 

the following system 

fj.(p ~ 2div D{(p) = F{ip) in S (44) 

^ = C ondS (45) 

admits a unique solution ip in L^(0, oo; H^(iS) H H^(0, oo; H^(5)), for ^ > large 
enough. This function if satisfies the conditions 



(f ■ ndT = 0, / ifdy = 0, / y A ipdy — 0, 

dS JS JS 

and there exists a positive constant C > such that 

ll'yi'llL2(0,oo:H3(5)nHi(0,oo;Hi(5)) < C*!! Cli L2 (0.oo;H5/2(5))nHi (0,oo;Hi/2(5)) • (46) 

Besides, if (1X2 G L'^{0,oo;li^/'^{S)) ll^O,oo;'H.^/'^{S)), then the solutions ipi 
and ip2 associated with Ci and C2 respectively satisfy 

\\f2 - <<5l||L2(0,oo;H3(5)nHi(0,oo;Hi(5)) < C'IIC2 - Cl II L2(0,oo;HV2 (5))nHi (0,oo;Hi/2(5)) • 

(47) 

For if e L2(0,oo;H3(5)nHi(0,oo;Hi(5)) obtained from C, we now define X^ - 
Ids e Wa(S'So) as follows 



Xli;t)^y+ e-^V(-,s)ds 4^ ^ ^ ^ 

X*{;0)=lds. 
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Thus the esthnates (46) and (47) become 

llX^-Id^ll < C||C||L2(0,oo;H5/2(5))nHi(0,oo;Hi/2(5)), (48) 
ll"'^C2 ^ ^Ci II - '^llCa - Cl|lL2(0,oo;H=/2(5))nHi(0,oo;HV2(5))- (49) 

The resuh of this proposition could be actually reduced to saying that there exists a 
linear continuous operator from L2(0, oo; H5/2(5)) nH^ (0, oo; li^^^{S)) to Wa(S'JL): 
but the means that we use for obtaining it (that is to say by considering system 
(44)-(45)) is actually useful for a key point in the second part of the proof. 

5.3. Second step of the proof. Let us consider a control G WxiS^) which 
has been obtained in the previous subsection from a boundary velocity ^. Instead 
of projecting on a set of controls satisfying the nonlinear constraints required by 
Definition 2.1, we prefer projecting the displacements X^ — Idg, because we choose 
the space yV\{S^) as an Hilbertian framework. We denote the displacement by 

Z*^ = X*^-Us. 

The goal of this subsection is to define (in a suitable way) a mapping X* S 'W\{S^) 
which satisfies the nonlinear constraints of Definition 2.1. We associate with it the 
displacement 

Z* = X* - Ids, 

so that the wanted mapping is now Z*. We can decompose such a mapping as 
follows 

Z* =z^ + {z* - zi) = zi + {X* - XI) . 

Let us define the differentiable mapping 

5: Wx{S%,) ^ Hi(0,oo;M3) X Hi(0,oo;R3) X Hi(0,oo;M) 
Z* ^ 'S{Z*) 

^iZ*){y,t)^( /(Z*+Id5)A— ; / (cof (VZ* + I^a))^ 



and the spaces 



dt ' Js dt ' Jqs dt 

Sni = {Z* eWx{Sl)\UZ*) = 0}, 
El = {Z* e Wa(5So) I ^o;?(^c) = 0} ' 



, , dZ* r dZ* r dZ* 



where 

Note that £i is a space where lie Id^ and X^. That is why the constraints satisfied 
by Z^ and X^ arc the same. Note that the space £„; takes into account the nonlinear 
constraints adapted to the displacements Z*. The purpose of this paragraph is to 
project any displacement Z^ € £i on the set £ni, provided that the displacement 
is close enough to 0. 
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Figure 1. Projection of a given mapping on the set of admissible displacements 



The definition of such a projection is given by: 

Theorem 5.2. Let be e Wx{S?^). If Z^ is small enough in Wx{S'^), then there 
exists a unique mapping Z* G £ni such that 

Moreover, we have that 

\\Z* - Z^\\w^(s°^) = o{\\Z^\\w^{s°^)) ■ (50) 

Thus we denote by V : Z^ Z* the projection so obtained. 

If the displacements Z^^ and Z^^ are close enough to in W\{S?^), then 

11(^2 - ^C2) ^ (^1 ^ ^Ci)IIwa(S5„) < ^ ^CiIIwa(SO,) 

xA-(||Z*J|wa(5S„) + II^cJIwa(5S„)): (51) 
with Zl = VZ'^^ and = VZ^^^, and K{r) — > when r goes to 0. 

Remark 6. Note that in this statement we do not need to assume that Z^? e 
But the way we have constructed X'^ such that Z^? = X'^ — Ids G in the previous 
subsection will be useful in the proof below. 

Proof. The proof of this theorem is an application of Theorem 3.33 of [3] (page 74), 
that wc state as follows : 

Theorem 5.3. Let W be a Hilbert space, TZ a Banach space, and g a mapping of 
class from W to U, such that 5"^({0}) ^ 0. Let be Z^ G W, and Z^ e 5~^({0}) 
such that Dx^g is surjective. Then there exists e > such that if \\Z^ — Zq\\)^ < e, 
then the following optimization problem under equality constraints 

min 

admits a unique solution Z* . Moreover, the mapping Z^ h- > Z* so obtained is . 
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In order to apply this theorem with 

W = Wa(S'^), 7^ = HJ(0,oo;M3) xHJ(0,oo;R3) xHJ(0,oo;M), 

the only nontrivial assumption to be verified is that the mapping Dq^ is surjective. 
For that, let us consider 

(a, 6, c) e Hj(0, oo; R^) x HJ(0, go; M^) x HJ(0, go; R). 

An antecedent Z* of this triplet (a, 6, c) can be obtained as 



where ip is the solution of the following system 

flip - 2div D{ip) = F{ip) + Fa,b in S, 



c 



on 95, 



for /i > large enough, with 

F{Z*){y, t) = ^ 2i^(Z*)ndr^ + (^I^' J^^ Ids A 2D{Z*)ndr^ A y. 

The previous study of system (44)-(45) (see section 7 of Part I) can be straight- 
forwardly adapted to get the existence of a solution in ip e L^(0, od; H^(5)) n 
H^(0, oo; H^(iS)) for such a system, and thus a displacement Z* e 'W\{S^). 
Since the projection V : i-^ Z* so obtained is C^, we can notice that its differ- 
ential at is the identity ld>V;^(so ), and thus a Taylor development shows that 

Z* = Z*^+0 + o{\\Z*^\\w,iso^)), 

\\Z* - = O (||Z^||yi;^(50^)) . 

For Z^_^ and Z^^ close to 0, the estimate (51) is obtained in considering a Taylor 
development around Z^_^ for the mapping V — Idyy^ (sSo ) • 

{Z; - ZD [Zl ZD = [Dz^^V Idw.tss.)] {zi ZD 

+o(ll^C2 " ^Cillw;>(ss„)) ■ 
Since Z* i— > Dz-V is continuous at 0, we have 

DziV - Idw;,(sS„) ^ when ||Z^J|w^(so^) goes to 0, 

and thus we obtain the announced estimate. □ 

Then, from the displacement — X'^ — Ids we can define a deformation X* as 
follows 

X* = V{X^ -lds)+lds. 
This deformation is admissible in the sense of Definition 2.1. 

The interest of such a decomposition (that is to say X* = X^ + {X* — X^) with X* 
given by Theorem 5.1) lies in the fact that the admissible control X* so decomposed 
will enable us to stabilize the nonhomogeneous linear part of system (l)-(9) thanks 
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to the term (see the previous section), whereas the residual term {X* — X'^ 
satisfies the property (50), which leads to 

11-'^* - Ids II Wa(s;j„) < C\\X^ -lds\\w^(si,), 



dt at 



= o(||X^* -Idsllw^fgo^)) 

L2(H5/2(as))nHi(Hi/2(as)) 



By combining the second and the third inequality to the estimate (48), we then 

obtain respectively (40) and (41). 

Lastly, the estimate (51) is reformulated as follows 

K {\\Xl - Idsl|w.(ss„) + II^C2 - Id5||>v.(5S.)) X 11^0 - ^cJlw.(5!L). (52) 

/ V / L2(0,oo;H5/2(a5))nHi(0,oo;Hi/2(05)) 

K (ll^Ci - Ifl^llwAlSSo) + ll-'^G ^ Ws||>v^(50^)) X ||X^^ - X^J|w^(50^), 

(53) 

where K{r) — > when r goes to 0, and where 

XI = 7'(X*^ -Id5) + Id5, 

x; = 7'(x^; -id5) + id5. 

The inequality (52) combined to (48) and (49) leads to the estimate (42), and the 
inequality (53) combined to (48) and (49) leads to the estimate (43). 

6. A fixed point metliod. System (l)-(9) is transformed into system (23)-(29). 
Before proving the stabilization to zero, with an arbitrary exponential decay rate 
A > 0, of system (l)-(9), let us first prove the stability of system (23)-(29) for all 
A > 0, for some well-chosen deformation X*. 

6.1. Back to tiie nonlinear system written in a cylindrical domain. In 

system (23)-(29), the mapping X* has to be admissible (in the sense of Definition 
2.1). It has to be chosen also in order to stabilize the linear part of this system. 
For that, we decompose formally the function e'^^^^f- on dS as follows 

Let us choose in system (23)-(29) the function C, in the following feedback form: 

C = }C\{u — G{u^h' ,Cj),h' ,iJj). 

Provided that C, is small enough, by Theorem 5.1 we can now define an admissible 
deformation X* . From this deformation we can define a change of variables X 
and the corresponding mappings X and Y (see (21)) which enables us to define 
the right-hand-sides of system (23)-(29); More precisely, we rewrite this system as 
follows 

du 

— \u-'iyAu + \/p = F{u,p,h',Ld), inj"x(0, c»), (54) 

dt 

div -u = div G(it, /i', w), inJ^x(0,oo), (55) 
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u = 0, onaOx(0,oo), (56) 

ii = h'{t) + u}{t) Ay + ICxiii - G{u, h' , w), h' , w) (57) 

+ (e^*^ - + W{u, h', {y, t) e dS x (0, cx)), (58) 

Mh"-\Mh' = - a{u,p)ndT + FMiu,P,h',Lj), in {0,oo), (59) 

JdS 

Iqlo' {t) — Xli^uj = — I y A a{u,p)ndT + Fj{u,p,u}), in(0, oo), (60) 

JdS 

{t(y,0) = uo(2/), 2/e /i'(0) - /ii e R^ w(0) = c^o e M^ (61) 

with 

F{u,p,h',Cj) = i^(L-A)M-e-^*M(M,^',w)-e~^*Nu-(G-V)p-e-^*tjAit, (62) 
G{u,h',Q) = {lRS~VY{X{y,t),t))u, (63) 

W{u,h',uj) = Co A {X* -Id), (64) 
C = lCx{u-G{u,h',Lo),h',uj), (65) 
FM{u,p,h\Lj) = -Me-^^ib Ah'it) 

-V j ^Vu(vf(X) -Irs) + (^Vf (X) -Irs^^Vu^^ vy(x)^ndr 

-/ ct(u,p) (vy(X) -lR3)^ndr, (66) 
Fi{u,p,uj) = -(r-/o)w' + A(r-/o)w-r w + e~^*rt2;Acj 

-i^y y A (^Vm (vf(l) - Irs) + (Vy(X) - lR3)^V{t^) Vf (l)^ndr 

yA (a{u,p)iVYiX) -lR3)nY dT 

OS ^ ' 



+ [ {X* - Id) A (a{u,p)\/Y{Xfn) dT. (67) 
JdS ^ ^ 



Note that the mapping X* is entirely defined by C, and so by the unknowns (u, h',oj), 
while the mappings X and Y are entirely defined by X* and the unknowns (/i',w), 
and thus by {u,h',ui). 

Remark 7. Note that the projection method used in order to define an admissible 
deformation X* from a boundary control ^ has been made in infinite time horizon, 
with regards to the functional spaces considered. It implies that the nonlinear 
system above is noncausal, that means that the control ^ chosen in a feedback 
form anticipate a priori at some time t the behavior of the unknowns for later 
times. In practice we could define a projection method for increasing times, but the 
corresponding Lipschitz estimates - in infinite time horizon - obtained above do not 
hold anymore. Anyway, let us show that this system admits a unique solution, and 
thus that it makes sense for all time. 



6.2. Statement. 
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Theorem 6.1. For (mo,/ii,wo) small enough in Hj^, system (54)-(61) admits a 
unique solution {u,p, h',uj) inR'^'\Ql^)xL'^{Q,oo;R\T))xR\0,oo;R^)xR^{Q,oo;M 
and there exists a positive constant C such that 

1I'"1Ih2.i(qo^) + ||p||l2(0,oo;H1(J^)) + II^'I1hi(0,oo;R3) + I|w|1h1(0,oo;R3) < C. 

6.3. Proof of Theorem 6.1. Let us set 

e = H2-i(g!^) X l2(0,oo;H1(J')) X h1(0,oo;M^) X h1(0,oo;R^). 

A solution of system (54)-(61) can be seen as a fixed point of the mapping 

N : H ^ H 

{v,q,k',w) 1-^ {u,p,h',u;) 

where h'^Lu) satisfies 
du 

— Au — lyAu + Vp = F{v, q, fc', ■uj) in x (0, oo), 

div ii ~ div G(u, fc', w) in J- x (0, cxj), 

u = Q on do X (0, oo), 

u = h'(t) +Lj{t) Ay + ICxiit- G{u,h',Lj),h',Ld) 



At 



e 



dX* 



- C] +W{v,k',w) on55x(0,oo), 



Mh"-XMh' = - a{u,p)ndr + FM{v,q,k',w) in(0,cx)), 
Jds 

/ow'(i) — A/qw = — / y A a{u,p)ndT + Fi{v, k' ,w) in(0, cx)), 

JdS 

u{y,0)^uo{y), y eJ^, h'{0) = hieM.\ w(0) - e 

with 

C — IC\{v ~ G{v,k' ,u7),k' ,w), 

and X* which is obtained from C by Theorem 5.1. The system above is actually the 
nonhomogeneous linear system introduced in section 4. In particular the estimate 
(39) gives 

II"I1h2-MQSo) + II^P|Il2(0,oo;L2(J^)) + 1|^'1|h1(0,oo;R3) + II II Rl (0.oo;R3) < 
Co (||wo||h1(.F) + \hl\R3 + |wo|r3 + ll^^lw,?, fc',^^7)||L2(o,oo;L2(.F)) 

+ \\G{v, k', ■n7)||H2.i(Qo^) + \\W{v, k', ^)\\L^(o,oo;m/-^dS))nii^o,oo-uu^{dS)) 



L2(0,oo;H3/2(05))nHi(0,oo;Hi/2(a5)) 

+ ||FM(w,g, A:',w)||l2(o,oo;R3) + \\Fi{v,k',m)\\i^2(^Q ,^.jg 



dX* 



(68) 
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6.3.1. Preliminary estimates. Some estimates given below are obtained by using a 
result stated in the Appendix B of [10] (Proposition B.l), and that we recall in 
Lemma 8.1. Let us also keep in mind the regularities provided by Proposition 4 for 
mappings X and Y(X). 

Lemma 6.2. There exists a positive constant C such that for all {u,p, h' ,Cj) in H 
we have 

II (A - L) u||^2(o,oo;L2(.F)) < C II"IIl2(o,oo;H2(.f)) >< (69) 
\\VY{X)VY{Xf - Ir3||loc(o,oo;H^(^)) + || Af(X)||Loc(o,oo;Hi(^))) , (70) 

||(V - G)p||^2(o,oo;L2(.F)) < C\\VY{X) -Ik3||l~(0,oo;H2(.F)) I|Vp||l2(o,oo;L2(.F)) ■ 

(71) 

Proof. The only delieate point that consists in verifying that /SY{X) in L°°(0, oo; H^(J^)). 
For that, let us consider the i-th component of /SY{X); We write 



Ay;(x(-,i),t) = trace (v2i;(l(.,t),t)) 



with 



y^Y,{X{;t),t) = (v(vK,(l(-,t),t))) Vf(X(.,t),0 

= (v (vy,(x(.,t),t) - Ir,.)) vY{x{-,t),t), 

and we apply Lemma 8.1 with s = 1, /i = and k = 1 to obtain 
||AY;(X(-,i),i)llHi(^) < C||Vy(l(-,t),t)-lM3|jH2(^)||Vf(X(.,0,0llH^(^)- 

□ 

Lemma 6.3. There exists a positive constant C such that for all {u,p, h',uj) in H 
we have 

I|M^IIl2(0,oo;L2(.F)) ^ C'II"IIl2(0,oo;H2(.F)) ( II VF(X) - Ir3 | j l- (o,oo;H2 (.F)) + 1 

dX 



1|/i'IIl°=(0,oo;R3) + II'^IIl°°(0,oo;R3)I|-'^IIl'=(0,oo;H1(.7=')) + 



dt 



L°°(0,oo;Hi(.F))y 



^"IIl2(o,oo;L2(J^)) ^ II"IIl'=(0,oo;H1(.?')) H "H (0,oo;H2 (J^)) ^ 



(72) 



\VY{X) - Ir3 ||l=o(o.oo;H2(^)) + 1 j , (73) 

11^ A 'u||i2(o,oo;L2(J^)) < C'||w||Hl(0,oo;R3)||w|lL2(o,oo;L2(jr))- (74) 

Proof. There is no particular difficulty for obtaining these estimates. □ 

Lemma 6.4. There exists a positive constant C such that for all it G H'^'^((3j^) we 
have 

II^(^)IIl2(0.oo;H2(JF)) < C'll*i||L2(0,oo;H2(jr))||Vy(X) - Ir3 || L~ (0,oo;H2 , 
I|G(")IIh1(0,oo;L2(.F-)) < C (||u||Hl(0,oo;L2(.F))||Vy(X) - Ir3 1 1 L-= (0,oo;H2 (JF)) 

+ ll^*l|L2(0,oo;H2(Jf))||VF(X) - Ir3||w1.=^(0,oo;L2(J=-))) • 
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Proof. The quantity \/Y{X) lies in L^(0, oo; H^(J')). We apply Lemma 8.1 with 
s = 2, /X = and k in order to get 

l|G(u)(-,i)||H2(^) < C||ti||H2(^)||Vy(X(.,t),t)-/||H^(^), 
I|G'(w)IIl2(0,oo;H2(^)) < C'll'^llL2(0,oo;H2(jr)) 1 1 VF (X) - Ir3 || Loo (o,oo;H2 (^)) ■ 

For proving the regularity H^(0, oo; L^(J^)), we first write 

(vm)-i..)f .(|(vm,-/ 

The quantity Vy(X) lies in n^{Q,oo\l7{T)) ^ Wi-°°(0, oo; L2(J')), so that we 
have the estimates 



dG{u) 



dt 



L2(0,oo;L2(jr)) 



< C ||VF(X)-Ik3|1loo(o,oo;L-(^)) 



du 
'dt 



L2(0,oo;L2(J=-)) 

+ ^ Ik3|Ih'1'~(0,oo;L2(J^))||m|!l2(0,oo;L°=(J^)) 



□ 



Lemma 6.5. There exists a positive constant C such that for all {u,p, h',uj) in . 
we have 



ll^('^)llHi(0,oo;H3/2(a5)) ^ C'II'^I1h1(0,oo;R3)||-'^* - Wjl y^;(QO^ ) , 

FM{u,p,h',u}) < 

L2(0,oo;M3) 

C (^||/i'||l2(0,oo;B3)|1'^||h1(0,oo;K3) + ( 1 1 "1 1 (0,oo;H2 (J^)) + | (0,oo;Hl (J^)) ) 

(||Vf(X)|U<.(a5,(0,oo)) (l + l|Vf (X) - lR3|U<.(a5,(o,oo))))) , 

||i^7(u,P,w)|lL2(o^^;R3) < 

C ((1 + A)||/* — /o||L^(o,oo;R3)|jw||Hi(0,oo;R3) 

+ ll^*'llL2(0,oo;K9)ll'^|jL==(0,cx3;K3) + 1 1 1 1 L°° (0,oo;K9) 11^ II L°° (0,oo;B3) II || L2 (0,oc 
+ (ll"l|L2(0,oo;H2(:r)) + ||p||l2(0,oo;H1{J^))) X 

(||Vf(l) 

(||Vf (1) - /|U-(a5x(0,oo)) + \\X*- Id|U 



(75) 



(76) 



(77) 



with 



||-^*'||l2(o,oo;B5>) < C'|1-^*||l~(0,cx);L2(5)) 

- Io\\L°^{0,oo-m^) < C'||-Y*||l°°(0,cx);L2(5)) 
l|-f*llL=°(0,oo;R9) < C'l|-Y*|jL°°(0,oo:L2(5))- 



dX* 



dt 



L2(0,oo;L2(5)) 

dX* 



dt 



L2(0,oo;L2(5)) 



Proof. There is no particular difficulty for proving the other two estimates, if we 
refer to the respective expressions of W, Fm and Fj given by (64), (66) and (67). □ 
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For some radius R > 0, let us define the ball 
Br = i^{u,p,h',uj) e H I M = on (90, and 

II"IIh2,1(Q!^) + ||p||l2(0,oo;H1(J^)) + il^'llHi(0,oo;M3) + II II HI (0,oo;B3) < 2i?Co| 

with 

R = {\\uo\\Hi{r) + \hi\R3 + \ojo\r3) , 

and where the constant Co appears in the estimate (68). Note that Br is a closed 
subset of H. 

First, for R small enough, we claim that for (-0, q, fc', w) € Br the function 

C }C\{v — G{v,k' ,w),k' ,w) 

is small enough in L2(0, oo; H5/2(55)) n Yi^{fd,oo]'H}/'^{dS)), because of the esti- 
mates of Lemma 6.4 and the continuity of the operator 1C\ (see Part I, section 5). 
So by Theorem 5.1 we can define the corresponding admissible deformation X* , the 
mappings X and Y which stem, and thus - in virtue of Remark 4 - we can define 
properly the mapping M in Br. 

6.3.2. Stability of the set Br by the mapping Af. Let be {v,q,k' ,w) G Br, for R 
small enough. In the estimates provided by the previous lemmas 6.2, 6.3, 6.4, 6.5, 
note that from the estimate (93), (94), (95) of Proposition 4 combined to (40) and 
(41) of Theorem 5.1 we can deduce 





||F({},9, A:',tx7)||L2(o,oo;L2(.F)) 


- o{ 


liO 


Ihi(.f) 4 




^ k^O Rs) 






||G(w,fc',tI7)|lH2.1(Q0^) 


= o{ 


Uq 


Ihi(.f) 4 


- \hl\R3 H 


" k^O R3) 


\\W{v,k\ 




L2(H-V2(as))nHi(Hi/2(a5)) 


- o( 


liO 


Hi(Jf) 4 




- k^O R^) 


dt 


-c 




= o{ 


Uq 




- \hl\R3 H 


- Wo R3) 




L2(H3/2(a5))nHi(Hi/2(a5)) 














\\Fm{v-, q, k', to)||l2(o,oo:R3) 


= o{ 


Uq 


Hi(.F) 4 


- \hi\s.3 H 


- k^o R3) 






||F/(l>,fc',tI7)||L2(o,oo;R3) 


= o{ 


Uq 


Hi(.F) 4 


- H 


" k^O Rs) 



In particular, for R small enough the mapping Af is well-defined, and so we can 
write the estimate (68) that we recall 



Af{v, q, fc', w) 



< 



Co (^I|Wo||h1(.F) + \hl\R3 + |Wo|r3 + ||-F(v,g, fc',n7)||L2(0,oo;L2(j:-)) 

+ \\G{v, k', n7)||H2,i(go^) + ||M^('0, k', ro)||L2(o,oo;H3/2(a5))nHi(0,oo;Hi/2(g5)) 



L2(0,oo;H3/2(e5))nHi(0,oo:Hi/2(a5)) 

+ \\Fm{v, q, k' ,m)\\-^2(^Q^^.TS.3) + ||-F/(w, fc', t£7)||L2(o^oo;R3) 
Thus we have 



Af{v,q,k',w) <Co{R + o{R)) 



(78) 
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This shows that for R ~ (||uo||hi(j^) + l^ik^ + l^i^olRa) smaU enough the bah Br is 
stable by the mapping Af. 

6.3.3. Lipschitz stability for the mapping J\f. Let (wi, qi, k[, wi) and {v2, 92, ^21 '^2) 
be in Br. We set 

and 

U ~ U2 — ill, p = P2 — Pl, — 00 = 002—001, 

V = V2 — vi, q = q2 ~ qi, k' = k'2 — k'l, w = W2 — w\ . 

We also denote 

Ci = ^\[vi- G[vi,k!-y,wi),k'^,wi), 

C2 = ICxih - G{v2,k2,'UJ2),k2,'UJ2), 

and next Xi and X2 the admissible deformation given by Theorem 5.1 which are 
defined by Ci and C2 respectively. 

For R small enough the quadruplet {u,p,h',uj) satisfies the system 
Q'li 

— Xu — lyAu + Vp = F, in J" X (0, 00), 

ot 

div u = div G, in x (0, 00), 

u = 0, on do X (0, 00), 

u = h'{t) + u{t) Ay + JCx{u~G, h' , oj) + W, {y, t) e dS x (0, 00), 

Mh"-XMh' = - a{u,p)ndr + Fm, in(0,cx)) 

JdS 

IqU}' (t) ~ XIquj = — y A cT{u,p)ndT + Fj, in(0,cx)) 

JdS 

u{y^ 0) = 0, in J", h'{0) = € M^ w(0) = e M^ 

with 

F = F{v2,q2-,k'2,W2) ~ F{vi,qi,k[,wi), 
G = G{v2,k2,W2) — G{vi,k[,vji), 
W = W{v2,k2^W2)~W{vi,k[,mi) 

Fm = ^4/(1)2,92,^2:^2) - ^A/(wi,gi,fci,'!^l), 

Fi = Fi{v2,q2,~k2,W2) - Fi{vi,qi,k[,wi). 

The deformations and X2 induce respectively the mappings [Xi, Yi) and {X2, Y2) 
which - partially - define the right-hand-sides above. 

The right- hand-sides F, G, W, Fm and Fi can be expressed as quantities which 
are multiplicative of the differences 

V, q, k', 00, {x;-x:), (X2-X1), 
(vf2(^2) - VYi(li)) , (Ay2(^2) - Ayi(Xi)) . 
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For instance, the nonhomogeneous divergence condition G can be written as 

G = (vf2(l2)-Vyi(li)){)2 + (VYi(Xi)-lK3){;. 

Then the estimates of Lemmas 6.2, 6.3, 6.4, 6.5 can be adapted for this right- 
hand-sides, so that the estimates (93), (94), (95) of Proposition 4 combined to the 
estimates (42) and (43) of Theorem 5.1 enable us to prove that for R smaU enough 
we have 



L2(0,oo;L2(jr)) 


— o( 


(w, g, k' , m) 


h) 




= o( 


(v, q, k\ -nj) 


h) 


L2(H3/2(as))nHi(Hi/2(a5)) 


= o( 


{v, q, fc', w) 


h) 


-F'm L2(0,oo;R3) 


— o( 


{v, g. A;', tzj) 


h) 


P^/ L2(0,oo;R3) 


= o( 


(v, g. A:', tzj) 


h) 



Then the estimate (78) can be apphed for the quadruplet (u,p, h' ,Cj): 

\\(u,p,h' ,Uj)\\u < Co (||J^||l2(0,oo;L2(.F)) + I|G|1h2,1(Q0^) 

+ ll^llL2(0,oo:H3/2(a5))nHi(0.oo;Hi/2(a5)) 
+ 1|-P'a/1|l2(0,oo;R3) + I|-P'/||l2(0,oo;R3)) • 

Then we have for R small enough 

\\{u,p,h',uj)\\M < Co X o{\\{v,q,k',w)\\M) < i|| (w, g, fc', to)||h, 

and thus the mapping A/" is a contraction in Bj^. It admits a unique fixed point, 
that is to say there exists a unique solution of system (54)-(61). The announced 
estimate is easily deduced from what precedes. 

7. Conclusion. After having proven Theorem 6.1 for the unknowns {u,p, h' ,ui) of 
system (54)-(61), let us recall the relations given by (19)-(20) and (22) 

e^*u{x, t) = R(t){t(y(a;, t),t), e^*p(.T, t) = p{Y{x, t), t), 
e^%'{t) ^ Il{t)h'{t), e^^Lo{t) ^ R(t)t2;(t), 

and the relation given by (21) 

X{y, t) = hit) + Ii{t)X{y, t), {y, t) € F y. (0, oo), 

Y{x, t) = y(R(t)^(x - h{t)),t), {x, t) e Qoo = U -F(i) X {t}. 

t>0 

Let us deduce the stabilization of the main system in the sense of Definition 2.2. 

Theorem 7.1. For (uo,/ii,wo) small enough in H^(J') x K'^ x K'^, system, (l)-(9) 
is stabilizable with an arbitrary exponential decay rate, in the sense of Definition 
2.2. 

Proof. First, from Remark 3, we can see that for ('u,p, /i', Co) given the quadruplet is 
determined in a unique way; In particular we can define the rotation R associated 
with the angular velocity lo as the solution of the following problem 

^(R) = §(e-^'RLi)R = e-^^-REiu) -..^i-. ~T 

dt ^ ' with s(u;j ^3 — wi 

R(i = 0) = Ijjs, \— tli2 
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Now by considering 

{e^'h') = e-^*R (w A h') + Kh" , ^ (e^'u) = Ku' , 

we can deduce from the estimate of Theorem 6.1 that 

l|e^*^'||Hl(0,oo;R3) < C", l|e'^*w|lHl(0,oo;B3) < 

and by considering 



du 

-RV7i(2/,OVy(X(y,t),t) X 



(^e-^*A' + e-^*a;AX(y,t) + e^*^j 

we can deduce from the estimate of Theorem 6.1 and the estimates of Proposition 
4 that, since we have detVX = 1, the following estimates hold 

lle^*u||L2(0,oo;L2(Jc-(t))) = ||u||L3(0,oo;L2(J=-(t))) < C*, 

f(e^*u) < C. 

L2(0,c3o;L2(J^(t))) 

Finally by considering the equalities 

e^*Vp(X(2/,i),i) = R(t)Vy(X(y,i),i)^Vp(y,t), 
e^*Vu(X(y,i),t) = R(t)VM(w,t)Vy(X(j/,i),t)R(t)^, 

and for the i-st component of vector u the following equality 

e^*V2u,;(X(y,0,t) = V(Vu,;(X(y,0,t))Vf(l(y,i),i)RW^ 
we can use the estimates of Lemma 4 to conclude the proof with the estimate 

l|e^*PllL2(0,oo:Hi(Jc-(t))) < C", ||e^*w||L2(0,oo;H2(J=-(t))) < G. 

□ 

8. Appendix A: The change of variables. Let us consider an admissible defor- 
mation X* E yVxiS^) - in the sense of Definition 2.1 - which satisfies, in particular, 
for all i > the following condition 

BX* 

coiVX*{y,tf—-{y,t)-ndr{y) = 0. (79) 

as 

dX* 

The regularity considered for the datum e^* in this section is 

n{S°^) = L2(0,c»;H3(5))nHi(0,c»;Hi(5)). 

The goal of this subsection is to extend to the whole domain O the mappings 
Xs{-,t) and Ys{-,t), initially defined respectively on S and S{t). The process we 
use is not the same as the one given in [18]. Instead of extending the Eulerian flow 
given by the deformation of the solid, we directly extend the deformation of the 
solid, because the difference in our case lies in the fact that the regularity of the 
Dirichlet data - written in Eulerian formulation on the time-dependent boundary 
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dS(t) - is limited. 

The goal is to construct a mapping X such that 



detVX = 1, 

X = Xs, 



in J" X (0, oo), 
on dS X (0, oo), 
on do X {0,oo). 



8.1. Preliminary results. Let us recall a result stated in the Appendix B of [10] 
(Proposition B.l), which treats of Sobolev regularities for products of functions, 
and that wc state as: 

Lemma 8.1. Let s, ^, and k inR. If f £ H^+'^(J") and g G H*+''(J"), then there 
exists a positive constant C such that 

II/.9||h^h-^) ^ C'II/I1h=+''(.f)I|5||h=+»(.f), 

(i) when s + /i + k > d/2, 

(a) with IJ.>0, K > 0, 2s + + K > 0, 

(Hi) except that s + fi + K>d/2 if equality holds somewhere in (ii). 

A consequence of this Lemma is the following result. 

Lemma 8.2. LetlC* be in W{Qlo)- Then 

cofVX* e L°°(0,oo;H2(J-))nWi'°°(0,oo;L2(J-)), (80) 

and, if X — Idjr is small enough in VV((9!^), there exists a positive constant C such 
that 



jcofVAT - lR3||L~(H2)nwi'°°(L2) < C\\VX - Ir3 |iL==(H2)nwi 



(81) 



Besides, if X^ — Idjr and X2 — Idjr are small enough in VV(Qj^), there exists a 
positive constant C such that 



IcofVXj -C0fVXi||L-(H= 



)nwi 



= (L2) < C'llVXj - VXi||Loo(H2)nwi-°°(L2)- 



(82) 



Proof. For proving (80), the case d = 2 is obvious. For the general case, let us 
show that the space L°°(0, cxd; H2(J")) n W^'°°(0, cxd; L2(J")) is stable by product. 
For that, let us consider two functions / and g which lie in this space. Applying 
Lemma 8.1 with s = 2 and /i = k = 0, we get 

II/5||l°=(0,oo;H2(J^)) < C'||/|iL°=(0,oo;H2(.F))ll5llL~(0,oo;H2(.F))- 

For the regularity in W^'°°(0, 00; L^(J^), we write 



djfg) 
dt 



df ^,dg 



Using the continuous embedding ( J^) ^ L°° ( J^) , we get 



d{fg) 



dt 



< C 



L°°(0,oo;L2(jr)) 



df 



dt 

dg_ 

dt 



L°°(0,oo;L2(J^)) 



I15|Il°=(0,oo;H2(J^)) 



l/ll 



L=(0,oo;L2(J=-)) 



L°'(0,oo;H2(J=')) ' 



28 SEBASTIEN COURT 

and thus the desired regularity. Thus the space L°°(0, oo; H2(J'))nWi'°°(0, oo; 

is an algebra. The estimate (81) is obtained by the differentiability of the mapping 

VX* M- cofVX* (see [1] for instance); More precisely, we have 



lK3-C0fVX = (^VX -lK3j -div(^X -Id^j +0(^||X -ld^ll^(QO^) 

SO that we get 

j|cofVX - lR3|lL=(H2)nWi'°°(L2) < Cj|VX - lR3||L°=(H2)nWi' = (L2)- 

The estimate (82) can be obtained by the mean-value theorem, so its proof is left 
to the reader. □ 



8.2. Extension of the Lagrangian mappings. Let us first extend the deforma- 
tion of the solid X* to the fluid domain J^, in a mapping that we have already 
denoted by X . 

Proposition 2. Let X* — Ids G Wa(5J^) he an admissible deformation, in the 
sense of Definition 2.1. Let us assume that X* — Ids is small enough in yV\{S'i^), 
that is to say that the function 



{y,t) ^ e 



xtdX* 



dt 



is small enough m L^(0, oo; H'^(5)) HH^ (0, oo; H^(5)). Then there exists a mapping 
X* e yV{Q"oo) satisfying 




in T X (0, oo), 
on dS X (0, oo), 
071 do X (0, oo). 



and such that 



|X*-Id^|| 



< C 



-.At 



dX* 



dt 



L2(H5/2(a5))nHi(Hi/2(a5)) 



(83) 



(84) 



for some positive constant C independent of X* . Besides, if X^ ^Ids and X 2 — Ids 
are two displacements small enough in W\{Sl^^), then the solutions X^ and X2 of 
problem (83), corresponding to X^ and X| as data respectively, satisfy 



1^2 -^ll 



< c 



,At 



dx; 

dt 



dX^ 



dt 



L2(Hi5/2(a5))nHi(Hi/2(a5)) 



(85) 



Proof. Given the initial datum X*{y,0) ~ y ioi y G S, let us consider the system 
(83) derived in time, as follows 



( cofVX' 



dVX 
dt 



dX dX* 



9K 
dX 

~dt 



dt 



in J" X (0, 00), 
on dS X (0, 00), 
on do X (0,00). 
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This system can be viewed as a modified nonlinear divergence problem, that we 
state as 



dX dX 



fix' 



dt 
dX 

~dt 



dt 



= 



in J" X (0, oo), 
on dS X (0, oo), 
on do X (0, oo), 



(86) 



X(.,0)=ld^, 



or 

dt 



(•,o) = o. 



with 



f{X*) = (lK3-C0fVX' 



dt 



If we search solutions to this system which are continuous in space, in using the 
Piola identity we can verify that the compatibility condition for this divergence 
system is nothing else than the equality (79). 

A solution of this system can be viewed as a fixed point of the mapping 

miQlo) ^ Wa(OSo) 



(87) 



where X2 satisfies the classical divergence problem 



di^ 



dX^ 
Tdt 



fix,) 



dX^ _ dX* 

dxl 

dt 



dt 



in J" X (0, 00), 
on dS X (0, 00), 
on do X (0, od). 



X2(.,0)=Id^, 



dX^ 
"dT 



(•,0) = 0. 



Indeed, let us first ver ify that for X - Idjr e Wa(QSo) wc have e^*f{X ) e 
L^(0, 00; H^(J^)) n 11^(0, 00; L^(J^)). For that, we recall from the previous lemma 
that cofVX* € L°°(0,oo;H2(J')) n W^'°°(0, c»; L2(J")), and wc first use the result 
of Lemma 8.1 with s = 2 and = k = to get 



l!e'^7(^ )IIl=(h^(^)) < C||lM3-cofVX |1l~(h2(^)) 

For the regularity in 11^(0, 00; L^(J^)), we write 

' .d^WX*\ dcoiWX* 



dVX 



dt 



L2(H2(.F)) 



.xt dfiX) 
xtdfiX*) 



cofVX 



.At 



df^ 



dt 



dVX 
dt 



dt 



< q|lK3 -C0fVX*|jL==(H 



(H2(.F)) 



L2(L2(^)) 



dVX 



dt 



Hi(L2(^)) 



+c 



dcofVX* 




xtdyx* 




dt 


L°°(L2(.F)) 


dt 


L2(H2(.F)) 
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where we have used the contmuous embeddmg H^(J^) ^ L°°(J^). Thus there exists 
a positive constant C such that 



< C|llR3-C0fVX |lL°°(H2(J^))nWi-°=(L2(J^)) 



(9i 



The estimate (88) shows in particular that the mapping T is well-defined. More- 
over, for the divergence problem (88) there exists a positive constant C (see [8] for 
instance^) such that 



dt 



< C ||e^*/(^i)llL^(H^(^)) + 



L2(H3(.F)) 







at 


L2(HV2(as))) 



and also 



dt 



< C (l + A)||e^7(^i)liHML^(^)) + 

Hi(Hi(.F)) V 

Thus there exists a positive constant Co such that 
xtdXl 







dt 


Hi(Hi/2(a5))y 



dt 

Let us consider the set 



< Co ||e'^7(^i)||H^,i(Qo^) + 







at 





(89) 



"Br = 



{Z* eWxiQ'L), \\Z*\\wm'L)^R} 



with 



R = 2Cn 



At 



dX* 



dt 



Notice that a mapping Xi —Idjr g *B/j satisfies in particular the following inequality, 
obtained in the same way we have proceeded to get the embedding (15): 



l^i-Id^l 



< C 



At 



dX^ 



dt 



= ciix;-id^i 



< CR. 



Then the inequality (89) combined to the estimates (88) and (81) show that for 
Xi — Idjr £ 05 _R we have 

ll^2-Id^llw.(Qo^) < Co(^Ci?2(Ci?+l) + J-^, 

and thus for R small enough, *Bi^ is stable by T. Notice that ^Bij is a closed subset 
of VV((5!^). Let us verify that T is a contraction in *8i^. 



■^Using results of [8], the nonhomogcncous Dirichlct condition can be lifted (see Theorem 3.4, 
Chapter II) and the resolution made by using Exercise 3.4 and Theorem 3.2 of Chapter III 
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For Xi ~ Idjr and X2 — Idjr in ^r, we denote Z 
which satisfies the divergence system 



TiXl - Id^) - %{X, - Id^) 



div 

dT 

dZ 
dt 



dT 
dt 



f{Xl)-f{X\) 



= 



in J" X (0, cxd), 
on dS X (0, 00), 
on do X (0,cx)). 



and thus the estimate 
,dZ 



dt 



For tackhng the Lipschitz property of the nonlinearity, we write 

dWXl 



cofVXa - cofVXi 



dt 



(1^3 — cofVX^ 



-*x . dV{X,-X,) 



dt 



By reconsidering the steps of the proof of the estimate (88) and by using (82), we 
can verify that for R small enough the mapping T is a contraction in *B/j. Thus T 
admits a unique fixed point in ^r. 

For the estimate (85), if X^ and X2 are two solutions corresponding to Xi and X| 
respectively, let us just write the system satisfied by the difference Z = X2 — X^^: 

ffy* 

div -^^ fixl)-fiXl) in ^x (0,00), 



_ dt 
dZ' dXi 



dZ 
dt 



dt 



dXl 
"dT 







on dS X (0, 00), 
on do X {0,00). 



Then the methods used above can be similarly applied to this system in order to 
deduce from it the announced result. □ 



Let us now consider h G H^(0, 00; M'^), and R e H^(0,oo;M^) which provides 
UJ G H^(0,oo;R'^). Let us construct a mapping X such that X{-,0) — Idjr and 

dctVX = 1 in J" X (0,00), 

X = h{t) + R(t)X* on dS X (0, 00), 
X ~ Idgo on dS x (0, 00). 

We cannot solve this problem as we have done for problem (83), because the proof 
would require the unknowns h and ui arbitrarily small enough, a thing that we 
cannot assume, even a posteriori. Instead of that, we utilize the mapping X 
provided by Proposition 2, and we search for a mapping X such that 

X{;t) = x''{;t)oX*{;t). 
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Such a mapping X has to satisfy 

detVX'^ = l inJ"x(0,oo), 
X = /i + RIda5 on 95 X (0,00), 
X^^Idgo onaOx(0,oo). 

For that, let us proceed as in [19]: We consider a cut-off function ^ G C°°(J^), such 
that ^ = 1 in a vicinity of dS and ^ = in a vicinity of dO. We define the function 

^R{x,t) ^ \h\t) A (a; - h{t)) -\\x- h{t)\''Lo{t), 

so that curl(5'fl)(a;, t) = h'{t)+uj{t)/\{x—h{t)), and we construct X^ as the solution 
of the following Cauchy problem 

-^{x*,t)=cnT\iC:SR)(x"ix*,t),t), X'^ix*,0)^x*, X* eX*{J^,t)^J^*. 

(90) 



We can verify (see [19] for instance) that the mapping X so obtained has the 
desired properties, and thus we can set 

Xiy.t) - x''(X*iy,t),t), (y,t) e J-x (0,00). (91) 

X (•,*) and X {■,t) are invertible, the mapping X{-,t) is invertible, and we 
denote by Y{-,t) its inverse. The mapping X presents the same type of regularity 
as the mapping X . We sum its properties in the following proposition. 



Proposition 3. Let X* be an admissible control - in the sense of Definition 2.1 - 
and X ttie extension of X* provided by Proposition 2 (for e^*^^ small enough in 
7^3(5*^^) ). Let X be the mapping given by (91). For all t > 0, the mapping X{-, t) is 
a -diffeomorphism from O onto O, from dS onto dS(t), and from T onto J^{t). 
We denote by Y{-,t) its inverse at some time t. We have 

{y,t)^X{y,t)eW{QL). 

detVX(j/,i) = 1, for all {y,t) e J" x (0,oo). 

The proof for the regularity of X can be straightforwardly deduced from Lemma 
9.3 in the Appendix B of this chapter. We do not give more detail in this section, 
because here the aim is only to get a change of variables which enables us rewrite 
the main system as an equivalent one written in fixed domains (see section 3). 

9. Appendix B: Proofs of estimates for the changes of variables. Let us 

recaU that for X* — Ids & Wa(5^) Proposition 2 enables us to define the extension 
X* e WxiQ"oo) + Id^ satisfying 

detVX* = 1 in J" X (0,00), 

X* ^ X* ona5x(0,c»), 

X* =ldao on 50 X (0,cx)), 

and 



dX* 



dt 
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For h G H2(0,oo;M3) and R e H2(0,oo;M9) which provides lo e H1(0,oo;R3) such 
that 



R(0) = Ik3, 













-UJ2 


UJl 







we can define X through the problem 

—-(x*,t)=cnTmR){x"ix*,t),t), X"ix*,0)=x*, X* &X*{F,t)=F*{t), 
at 

where ^ is a regular cut-off function, and 

dii{x,t) = ]^h\t) A (a; - h{t)) -\\x- h{t)fio{t). 



Then we define 



and 



— R — * 

X = X oX , 



X = IL^{X - h). 



Lemma 9.1. For X^ — Idg and X| — Ids small enough in Wa(S'J^), let X\ and 
X2 be the solutions of Problem (83) (see Proposition 2) corresponding to the data 
X* and X2 respectively. If we denote by Y^^-^t) and Y2{-,t) the inverses of Xi 
and X2 respectively, we have 



^^2(^2) - Vyi(Xi)||L°=(0.oo;H2(:r))nWi.°°(0,oo;L2(jr)) 



at 



dt 



(92) 



Proof. Let us recall the estimate (85), obtained for X^ — Ids and X| — Ids small 
enough in >Va(S'!^): 



< C 



dt 



\\^2 ^illvv(Q!!,) 
First, let us give two intermediate estimates: 

- lR3||Loo(H2)nWl'°°(L2) < 

and 

11^5^2(^2) - lR3|lL°=(H2)nWi'°°(L2) < 



At 9X2 xt 



At 



dt 



dX* 



dt 



^(■sSo) 



L~(H2)nWi.°°(L2) 



1 - CIIVX2 - lR3||L°°(H2)nWi.°°(L2) 

< 2IIVX2 - lR3|jL=c(H2)nWl'°°(L2), 

provided by the equality 

vf;(x;) - = (ik3 - vxi) (vf;(x;) - 1^3) + (i^a - vx;) 
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and by the fact that L°°(0, oo; H2(J")) n Wi'°°(0, oo; L2(J')) is an algebra. Then 
the estimate (92) is obtained by writing 

(i + iivf;(x;)-im3||) iivx;-vx;ii 

1-C||VX;-Ik3|| ■ 

□ 



yYl{xl)-yY\{x\) = 



\\VYI{XI)-VYI{XI)\\ < 



Lemma 9.2. Let Xi and X2 be the extensions defined by problem (90), with data 
(/ii,Ri) e H2(0,cx);M3) X H2(0,oo;R9) and (/i2,R2) S H2(0, c»; M^) x H2(0, 00; R^) 
respectively. Then we have 

\\^2 ~ -'^1 IIh2(0,oo;H'*1-°=(R3)) < O (^||/l2 - /1i||h1(0,oo:R3) + ll'^2 - '^i1|h1(0,oo;B3)) > 

where we recall that 

h[ = e^*Rlh[, h'^ = e^*R^/i^, 



Proof. The change of variables given by a mapping X is slightly the same as the 
one utilized in [19]; In considering the writing 

R{tfdR{x,t) ^ h'{t) A Il{tf{x - h{t)) - - h{t)\'6j{t), 

the steps of the proofs of Lemmas 6.11 and 6.12 of [19] can be then repeated, with 
the difference that in infinite time horizon we rather have 

\\^2 - -'^1 I|h2(0,oo;W.°°(R3)) < J'^''^' (^||/i2 - ^ii|Hl(0,oo;R3) + !|'^2 - <^i||h1(0,oo;R3) 

+ ||/l2 - /lli|L==(0,oo:R3) + 111^2 - Ri||l''=(0,oo;R9)) , 

where is bounded when hi, /12 are close to and Ri, R2 are close to Ijjs. In 
order to estimate ||/i2 — ^1 |1l°°(o,oo;R-^) and j|R2 — Ri!|l°°(o,oo;R3)j '^e first apply the 
Gronwall's lemma on 

-|(R2-Ri) = (R2-Ri)§(c:'2) + RiS(c:'2-(:^i) 

ot 

(R2-Ri)(0) = 
in order to get 

||R-2 - Ri!|l°=(0,oo;R9) < C'||£D2 - tJl || Ll (0,oo;R3) exp (C|| (D2 || Ll (0,oo;R3) ) • 

Besides, it is easy to see that 
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SO that IIR2 - Ri||l~(o,oo;R9) is controlled by \\u)2 - '^i||l2(o,< 
11^2 — ^1 |iL=°(o,oo;R3) can be treated in writing 

||'l2 — /il||L=°(0,oo:K3) < 1 1 '^-2 ^ ^1 1 1 Li (0,oo;B3) > 

= (R2-Rl)/l'2+Rl f/lj,-/l0: 



Then the term 



h'^ - h 



11^2 ~ '4IIl1(0,oo:1 
II^2||l1(0 

11/7,2 — /li||Li(0,oo:S 



< j|R2 - Ri||l=(0,oo;R3)II'i2|Il1(0,oo:R3) + 11^2 ~ ^'l II Li (0,oo;R3) 

< ^=II/i2||l2(0,oo;R3), 



< 



2A 
1 



11^2 ~ /ill|L2(0,c 



Finally, it is easy to verify that 

Il'i2 ^ 'illlHl(0,oo;R3) < (1 + A)||/l2 - /i'i||h1(0,oo;R3)j 
11(^2 - '^i||h1(0,oo;R3) < (1 + '^)||'^2 - '^i||h1(0,oo;R3)- 



□ 



Lemma 9.3. Let Xi and X2 be defined by 



x, = xfoxl, 



X, 



— R — * 

X^ 0X2, 



* * R R 

where X-^, X^, X^ and X2 are given in the assumptions of the previous lemmas. 
Then we have 



1^2 -Xi! 



rK{r), 



where 



At 5X2* 



dXl 



dt 



dt 



^(■sSo) 



^ — 11/^2 - /lillHi(0,oo;R3) + ||'^2 - ^^iUhMO^ooiRS) + 

and K{r) is bounded when r goes to 0. 
Proof. Let us write 

— R — * — R — * / — R — R\ — * 

X2 — Xi = X2 o X2 — X2 O X J + I X2 ^ X-^ \ o X^. 

— R — * — R — * 

For tackling the difference X2 ° X2 — X2 ° Xi, let us apply Lemma A. 3 of the 
Appendix of [4]; We get the estimate 

l|Ar2 o a:2 - ^2 o a'i||h3(j:-) < c\\X2\\ci(T) 



^2 - -'^iIIh3(j^) X 

(||^;IIh3(^) + II^;||^3(^) + i), 



and thus the regularity in L°°(0, oo; H3(J")). The regularity in Wi'°°(0, oo\ Hi(J")) 



can be also obtained in applying Lemma A. 3 of [4] for the time derivative of X 
^2 ~ ^2 ° ^1- 

For the term {X2 — -'^i I o X^^ wc apply Lemma A. 2 of the Appendix of [4]; We 
get the estimate 



< 



I-'^iIIh-' 



1 



H3(^) 

and thus the regularity in L°°(0, 00; H'^(J^)). Here again the regularity in the space 
W^'°°(0, 00; H^(J')) is obtained in applying the same lemma on the time derivative. 

□ 
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Proposition 4. Let Xi and X2 be defined by 

Xi = nJiXi - hi), X2 = R^(Xi - h2), 

where Xi and X2 are given in the assumptions of the previous lemma. Then 

ll^2-^ilU(Q0^) < rK{r), (93) 
||Vy2(^2)- Vyi(Xi)||Loc,(H2)nwi.=o(L2) < rK{r), (94) 
||Ay2(X2)- Ayi(Xi)||L»(Hi) < rKir), (95) 

where 

f = ||/i2 ~ '*il|Hl(0,oo;R3) + ||li2 - 'il||Hl(0,oo;R3) + 

and K{r) is bounded when r goes to 0. 
Proof. For proving (93), it is sufficient to write 

X2-X1 = R|^(X2-Xi) + (R^-R[)(Xi-;ii)-Rj(/i2-/ii) 
and to apply the previous lemma. The estimate (94) can be proven exactly like 
the estimate (92). Finally, for the estimate (95) we denote by Yi^i and 1^,2 the i-th 
component of Yi and Y2 respectively, and we write the equality 

v2y;,2(^2) - v2f,,i(Xi) = (v (vy;,2(X2) - vfol^i))) vf2(^2) 

+ (v (y,,i(Xi))) (vy2(^2) - vfi(Xi)) , 

and apply Lemma 8.1. □ 
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